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We argue that the occurrence of late-time acceleration can conveniently be described by first- 
order general relativity covariantly coupled to fermions. Dark energy arises as a Bardeen-Cooper- 
Schrieffer condensate of fermions which forms in the early universe. At late times, the gap and 
chemical potential evolve to have an equation of state with effective negative pressure, thus naturally 
leading to acceleration. 
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I. INTRODUCTION 

Current cosmological observations point to a universe dominated by a negative-pressure fluid component, dubbed 
dark energy, whose origin is unknown. If the equation of state of this fluid does not evolve and is w = p/ p = 
then dark energy is a cosmological constant. In this case, observations give a value of the cosmological constant 
that is 120 orders of magnitude smaller than the theoretically expected evaluation. We still lack a convincing way 
of understanding this issue. As a result, an approach towards making progress is to assume that the cosmological 
constant/dark energy stem from some new physics0 

Ever since the observational evidence of late-time acceleration, model builders have sought to find a candidate for 
dark energy. However, this task is daunting as it is hard to identify an existing degree of freedom in the standard 
model or general relativity that (i) has negative pressure, (ii) is homogeneous on horizon scales and (iii) matches the 
observed energy scale. This initially led to quintessence models where a new fundamental scalar degree of freedom 
with a fine-tuned potential can be adjusted to yield a late-time acceleration tracking dark matter. Alternatively, 
investigators have invoked infrared modifications of general relativity, for instance f{R) or Gauss-Bonnet gravity. 
While successful, many representatives in both classes of models suffer from fine-tuning or other theoretical problems 
due to the introduction of new degrees of freedom. 

In this work, we take a minimalistic approach to dark energy assuming no extra degrees of freedom except fermionic 
matter on a flat Friedmann-Lemaitre-Robertson-Walker (FLRW) background. A finite density of fermions in the 
early universe can undergo a Bardeen-Cooper-Schrieffer (BCS) condensation due to a covariant attractive channel 
from general relativity. The system is described by a set of transcendental equations that give relations between the 
scale factor, the fermion gap and the chemical potential. It was already shown in j^-HI that such a BCS condensate 
can play an important role in the early universe by resolving the big bang singularity via a bounce. Remarkably, we 
find that depending upon some of the parameters of the theory, the same condensate can also affect the history of 
the universe at late times. In this paper we will show that the nonperturbative potential of the fermion gap can lead 
to late-time acceleration^ 

This is achieved by extending the analysis of 6], which was applied near a cosmological bounce where the gap 
equation could be obtained in a Minkowski spacetime. Here we work on a FLRW background and consider the 
evolution of the universe from the bounce on. As a consistency check, we obtain gap equations similaiH to those 
in and reproduce the same bouncing cosmology. A mechanism of fermion condensation was argued to be relevant 
to dark energy also in 0, @ (see also [3, [l^ ) ; for another approach describing fermions and condensates on FLRW 
spacetimes, see [TTI - [l3| . 

The paper is organized as follows. In Sec.|ll]we introduce the BCS mechanism on a FLRW background and derive 
the effective equations of motion. This section is a little technical and the reader mostly interested in the cosmology 
can skip ahead to Sec. IIIIl where we classify the cosmological solutions and find numerical examples which accelerate 
at late times. Section |IV] is devoted to discussion. 



II. COSMOLOGICAL BCS THEORY 

We use Greek indices /i, i', • • • = 0, . . . , 3 for spacetime directions in a nondegenerate manifold M (with signature 
— hH — h) and capital Latin indices I, J, ■ ■ ■ 0, . . . , 3 for the internal Lorentzian tangent space T^. Spatial directions 
on M and Tj^ will be denoted as a,b^c, . . . and i, j, fc, . . . , respectively. Dotted and undotted indices a, a label spinor 
components. We work in units h = c — G — 1. 

We begin by studying the cosmology of first-order general relativity covariantly coupled to Dirac fermions. In such a 
system a four-fermion interaction emerges when we solve for the torsion. In what follows, we show that this interaction 
in a FLRW background realizes a BCS condensate whose potential generically leads to late-time acceleration. 



^ For alternative approaches which try to avoid dark energy by invoking large scale inhomogeneities see, for instance, [JQ. 

^ On the other hand, the present proposal does not address either the smallness or coincidence problem for the cosmological constant, 

although it does relax the former. These issues will require a better understanding of the regularization mechanism (see below), which 

goes beyond the scope of this work. 
^ Some of the coefficients differ slightly. This is mainly because here for simplicity we only consider the scalar channel of the four-fcrmion 

interaction, while in Q we also included the pseudoscalar interaction. Also, unlike in Q, we do not throw away the contributions from 

the antifermions. These differences do not affect the physics, since the coefficients have the same sign. 
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A. Fermions in FLRW background 

Let us start by considering pure general relativity described by the Hoist action Afterwards we sketch how 
the attractive four-fermion interaction emerges after solving for the torsion. The action on a curved manifold A4 is 

where is the gravitational field (vielbein, tetrad), e = | det ej\, -R^/'^ is the curvature of the spin connection A^'^ , 

7 is the Barbero-Immirzi parameter and R^^'^ — (-^^kl ^^iv'^ /"^ is the dual field strength. 

The first term in Eq. pi.l[) yields the Palatini formulation of the Einstein-Hilbert action, the latter emerging when 
inserting the solution to the associated equation of motion 



I. J - J [ 



(where F^^ are the Christoffel symbols) into Eq. and using g^ii, = e^rjij. The second term is identically zero 

on half-shell, due to the Bianchi identity for the Riemann tensor. OJ^^ is the metric-compatible, torsion-free spin 
connection. 

The story changes when we covariantly couple chiral fermions to the gravitational action. The Dirac action for 
fermions is 

= / d^xe (V;7^e/^V^V + cc.) , (II.3) 
where cc. denotes the complex conjugate and the covariant derivative is defined as V,^ = dp, — Clp, where 

= ^Aijpj:" , (II.4) 

and 

(II.5) 

are the Lorentz generators for spinors. 7^ are the usual Dirac matrices in Weyl (chiral) basis, 

where = (I, — cr*), = (I, cr') and cr* are the Pauli matrices: 

i_/Ol\ 2_ f -i\ 3 /lO 



^-[10)^ ^-Uoj' ^no-ij- ^"-'^ 

Note that a tetrad-based formalism is essential for the inclusion of fermions in the theory, since Dirac spinors live 
naturally in SU{2). 

The equation of motion for the total action S'h-I-5'd are solved in terms of a connection Aj'^ having two contributions 
the torsion- free spin connection and a torsion term related to the axial fermion current: 

V = <^N + c/^ (11.8) 

where Cj'^ is the tetrad projection of the contortion tensor, 

C,'-' = C;''elei] . (II.9) 

Square brackets denote antisynimetrization, X^^jj] ~ {Xjj — Xjj)/2. 

On solving for C^^"^ in terms of the fermionic field and inserting the resulting expression for Aj'^ in the total action, 
one obtains the four-fermion interaction (see the Appendix for a detailed derivation) 

Q /■ ^4 J5lJi 1 371 7^ 
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where Jg = 'ipl^l'^i' is the axial current. Therefore, the fermionic action is now 



'5'fcr — S-Q + S'int , 

where the covariant derivative in S'd is in terms of lo^'-' . For a flat FLRW Une element 



the vielbein is given by 



ds^ = -dt^ + a^{t) dxadx'^ , 



1 

a{t)6l 



where a{t) is the scale factor. The only nonzero structme functions and spin connections ujjjk are 
where H = d/a is the Hubble parameter. Then, 



aH 
Y 





-a^ 



(11.11) 



(11.12) 



(n.i3) 



(11.14) 



(11.15) 



Above we have assumed the gravitational action to be the usual Hoist action. However, when torsion is present (for 
instance, when it is generated by fermions or a spacetime-dependent Barbero-Immirzi field [T8l - [2^ ). it is natural to 
include it explicitly in the fundamental action, so that the Hoist term is completed by a torsion-torsion piece to form 
the Nieh~Yan invariant p7l[23j. 

There are several other reasons why to prefer the latter alternative. A second motivation is that the Hoist term is 
not topological and vanishes on half-shell, while one would expect to define the theory with topological contributions. 
Third, although mathematically correct the Hoist derivation does not respect the usual decomposition of torsion into 
its Lorentz irreducible components, and its trace part (a polar internal vector) turns out to be proportional to the 
axial current ^ [see Eq. (|A.8|) ]. 

In the Nieh-Yan case, the coupling A/^ = 2/(37r) no longer depends on the Barbero-Immirzi parameter. In either 
case the "bare" coupling is M > 1. As we will regard it as part of an effective coupling, our results will not be sensitive 
to the form of the classically vanishing part of the action. 



B. Weyl decomposition 



We are now ready to quantize the fermions on a FLRW spacetime. Since 



1 



the Dirac Lagrangian is 



(H.16) 



Just as in the case of Minkowski spacetime, it is convenient to decompose the fermions into two-component Weyl 
spinors. 



so that 



x'x) 



(11.17) 



(H.18) 
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One can also write the action in terms of left-handed Weyl spinors but there is now an extra term coming from the 
integration by parts (integration domain omitted from now on): 

since the x's are anticommuting Grassmann fields, x^X — ~XX^ ■ Foi' spatial derivatives, 

y" d^x o?x^a'd,x = - J d'^x a^d^x^'x'' = J d'^x a^x^'d^x^ , 

and Eq. (|II.18|) yields 



(11.19) 



where C = X^- The action is completely symmetric with respect to the particle and antiparticle Weyl spinors, ^ C. 
In other words, the expansion of the universe does not distinguish between particles and antiparticles. 

Moreover, since FLRW is conformally flat and the action is first-order in time derivatives, the latter reduces to the 
Minkowski action. To see this, let us first perform the conformal rescaling 



3/2, 



Then, 



p,UJ ■ 



Next we introduce the Fourier transforms 

f(x) = J d^pduj e-^["P-^-"*lCp,<^ , Cix) = J d^pdu; g-^I^P'^-^'lC, 
In terms of the Fourier components, the action becomes 

^g-i[ax-(p-p')+t(i^-i^')] 
d^pduj' dt ( — ^ Ue^^^,ep,c. + wCp,^'Cp,c. - '^l^'^"P^^p.uJ - Cp,^'^>.:Cp,c. 



-it(L0~L0') 



= (2^)3 / d'pdu^'dt k;^^,Cp,. +a.c^_^,Cp.. - 4,-''^>^?p." - CL'*>*Cp.. 



xe 



p 

it{uj — uj') 



Thus, the Dirac action in momentum space reads 

5d = (2^)^ I dS [c.(4^^ep,c. + C^l^Cp^u) - {^l^a^P^^p.^ + Cl.:^'p^Cp,u)] 



(11.20) 



(11.21) 



(11.22) 



(11.23) 



C. BCS condensation 

A simple and physically transparent way to understand the condensation mechanism is to introduce auxiliary scalar 
(gap) fields, which are proportional to the fermionic bilinears. The gap equation is then derived by integrating out 
the fundamental fermionic degrees of freedom. Our starting point is the four-fermion interaction term 



•Sint = / d'^xe 



•hiJl 



M2 



(11.24) 
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This term can be decomposed into a scalar, pseudoscalar and vector interactions using the Fierz identity 



(11.25) 



The last term is the higher-energy p-wave channel and we, as such, are going to ignore it. For simplicity, we will also 
drop the pseudoscalar condensate and only focus on the scalar one. Thus, our interaction reduces to 



S\nt — 



M2 



= d 



5'tlCC ) 



(11.26) 



where in the second equality we have introduced the auxiliary scalar A, which acts like a mass term for the fermions. 
For an FLRW background (spinorial indices restored). 



Sn 



(11.27) 



It is clear that a nonzero value for the auxiliary field A ~ -ipip would signal a (cosmological) BCS-like condensation. 
In order to find such a nontrivial value for A, one can take recourse to a mean- field approximation where the gap A 
is treated as a constant. With the same procedure of the last subsection, the mass term in momentum space is 



Combining it with the kinetic term (III.23p . 

mass 



(11.28) 



where ^ — S,a and ^ = The above can be written in four-component notation as 



-p, — 



p,—uj 



— p.— a; 



where Ap is a 4 x 4 matrix given by 



Ap — 



uj — a^pi A 
A Lu — a^pi 



(11.29) 



(11.30) 



(11.31) 



At this point we introduce a chemical potential fi in the action, which corresponds to having a nonzero number 
density of fermions. The matrix Ap is now modified to 



UJ — a'pi + fi A 

A u — a^pi — /i 



(11.32) 



The condition ^ < corresponds to a Bose-Einstein condensation of composite bosons. 



D. Effective action 



The resulting quantum theory is encoded into the path integral 

[PA][2?e][2?C]e*'^''='+^'"='=^ = / [PAJe'^"" « e^^-« U , (11.33) 



ISP ■ 

where we have integrated the Grassmann fields, defined the effective action ScB (often referred to as F in quantum 
field theory literature) and approximated the functional integral by the saddle point (mean- field approximation [24|] ) . 
The effective action Scs can be evaluated by performing the Gaussian integrals in terms of the fermionic coordinates. 
As usual, one ends up with a fermionic determinant. Eventually we have (see, e.g., (23. [25j) 

/d^ 
^^ln(dety4p). (11.34) 
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The determinant of Ap can be straightforwardly computed: 

det Ap = _ (|p| + _ ^2] [^2 _ (|p| _ ^)2 _ ^2] ^ 



(11.35) 



where A is the auxiliary field at the saddle point. This expression is not Lorentz covariant, as expected (see below). 
Accordingly, we are left computing 



eff 



r T 



(11.36) 



where 



/ = 



J 



(27r)4 



{ln[a;2 - (|p| + /x)' - A^] + ln[u;' - (|p| - - A^]} 



(27r)3 

= /l+/2. 



V(|p|+m)' + A2 + V(|p|-m)' + A2 



(11.37) 



The integral 7i is 



= ^ {p- i^rv^+A^ - £ ^ {p- i^rv^+A^ . 

To get I2, one has to simply replace /x O — /x in Ji. The second integral in Ii cancels the second in I2 and we are left 
with 



7= {p'+„^)^^T^. 

Jo 



The above integral can be regulated using the following formula: 



dp 



P 



r(i±^)r(i?-i±^) 



(p2 + ^2)B 2A2S-A-lr(B) 

Let us choose B = — 1/2 + £. Then we have {A = 2 and A = 0) 

1 



(11.38) 



(11.39) 



27r2 



r(3/2)r(£-2) ^ ^,^T{l/2)V{e-l) 



_A2(e-2)r(£- 1/2) A2(^-i)r(£- 1/2) 

xAFA2 



iA2r(£-2) + M'r(e-l) 



27r2r(e- 1/2)A2^ 

We expand up to 0{e) to obtain the effective action. Using expansion formulae such as 

1^ 



Tie 



-2^7? + 2 (7-2 + 2 ln2) V7re + 0(e^): 



(11.40) 



where 7 « 0.5772 is the Euler-Mascheroni constant, we have 



rA2 /I 



47r2[l - £(7-2 + 21n2)](l + £lnA2) 



A^ 
■47r2 



1 1 

e " 2 



21n2-lnA2 



- l + 21n2-lnA' 



Thus, the effective potential is given by 



M2^2 A2 

4 4-772 



— (-- i +21n2-lnA^ ] - ii^" [ - - l + 21n2-lnA 



4 Ve 2 



(11.41) 



(11.42) 
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In a renormalizable theory, the 1/e divergence can be absorbed using renormahzation conditions [2y, [23]- The four- 
fermion interaction term is nonrenormahzable in Minkowski and therefore the divergence cannot be ehminated. A 
standard approach is to interpret the regularization parameter in terms of a physical cutoff scale A, 



e 



(11.43) 



such that A remains finite. Here we just take a phenomenological approach and encode this arbitrariness, intrinsic to 
the model, in a reparametrization of the form A = e~^/^/2, where N is an 0(1) ^ O(IO^) free parameter: 



-A^-21n2. 



(11.44) 



The relation between e and A is really a matter of choice, so any sign and value of N is possible. If > then e < 0, 
which may happen if spacetime has a fractal structure in the ultraviolet. Later we will need N > 0(10'^) in order for 
the condensate to fit observations. Interestingly enough, such a range of values — 1 ^ e < may be compatible with 
a fractal interpretation of N, where the (early) universe shows an effective dimension slightly smaller than 4. Until 
we achieve a better control of the quantum theory, the issue of the physical interpretation of N will remain open, 
although we have just argued that it admits at least one possible accommodation. 
To summarize, the effective potential is 



M2 



4^ 



A2 

T 



N ■ 



In A" 



/i^ (iV + l + lnA^) 



(11.45) 



As already mentioned, the bare mass AI is 0(1) or larger. However, Eq. ([11.451) can be obtained also via a different 
regularization procedure which renormalizes the couplings of the theory [a, |24| | . This suggests that also M can be 
treated as a free parameter. It shall be our attitude in what follows. 

It is clear that the above potential has a minimum at dVcs/dA = given by the gap equation 



1 



= — \A'^(N + 1) - 2^^(N + 2) + (A^ - 2^^) A^l 



(11.46) 



If M is constant, this equation univocally specifies A as a spacetime function. It is useful to check that we recover 
the usual behaviour of the gap in the weak-coupling BCS limit [13, HI], where the fermion gas is diluted. For A ^ fi, 
Eq. pL46)) tells us that 



exp 



which is the familiar exponential suppression of the gap. 
The potential at the minimum is 



(H.47) 



a r 



167r2 



N - 



InA^ - Ai? 



(11.48) 



However, the potential energy that we have calculated includes the contribution from the chemical potential as well. 
The total number tiq of fermions is (23 



— d X ei/'7 ip — 



SSfc 



,AV 



'A^ + 1 + In A 



2^ -a^n. 



6fi dfi 2t:^ 

Assuming the system lies at the minimum of the potential, the total gap energy density of the fluid is given by 



(H.49) 



Pgap = Vinin + 

A2 - 



167r2 
A2 
327r2 



A^- 



InA^ I -4/x2(2A + 3 + 21nA2) 



[A'^ - 8fi^) {2N + 3 + 2 In A2) 



(H.50) 



We conclude this section by making two remarks, the first on the chemical potential. As soon as we fixed an FLRW 
background, we have chosen a "privileged" frame whereon one can define a homogeneous number density and a 
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chemical potential. Whenever we refer to these concepts, it is always with respect to this special FLRW homogeneous 
time slicing. Solutions of the field equations are free to break Lorentz invariance. One can have a chemical potential 
in a Lorentz- invariant Lagrangian, i.e., a term of the form ^Jivj'^ , while the background /^oi° is aligned in the cosmic 
rest frame. An example is the Kaplan-Nelson model of spontaneous baryogenesis. This is what happens also in big 
bang nucleosynthesis calculations, where one has to include the chemical potential of the different particle speciesQ 
This has nothing to do with breaking Lorentz invariance at fundamental level. Here we are assuming, in a self- 
consistent manner, that the theory and its false (perturbative) vacuum are indeed Lorentz-invariant, while the true 
vacuum is not when fJ. ^ 0- We will see that the measurable effect is in the cosmic expansion, and the later one looks 
into the evolution of the universe, the fewer the Cooper pairs one can detect, thus restoring relativistic physics at late 
times. 

The second comment is the following. We have introduced a physical cutoff by hand as in Minkowski nonrenor- 
malizable theories, but in [l^ it was argued that this type of BCS models in curved space may be renormalizable. In 
the same work, the running of the mass coupling M and a cosmological constant term was considered for a de Sitter 
background. Here we shall not endeavour to study the renormalization group flow for the particular, self-consistent 
cosmological background we will find later. However, related to the renormalization issue there is another. Although 
the tree-level interaction between fermions is attractive, the gravitational interaction could affect renormalization. 
This is in analogy with the phonon interaction in the standard theory, where the sign of the effective coupling gets 
flipped and one does not generically end up with a condensate [1^ [23] . There, one can consider the beta function 
for the running of the coupling M, and realize that the second-order (four-fermion) phonon-induced interaction is 
attractive (negative coupling) [28|. In our case, the graviton and the torsion are independent degrees of freedom and 
they can be dealt with separately. An attractive fermion coupling is generated by torsion after integrating it out. We 
may ask about torsion fluctuations, but they do not take place because torsion is not dynamical. Therefore, so long 
as we have an attractive channel at tree level and a Fermi surface, we will get a nontrivial gap (provided that the 
gap equation has a solution, which it does; see below). As far as torsion is concerned, the attractive nature of the 
tree-level interaction is guaranteed, unlike the phonon-induced screening (or giuon exchange) in typical scenarios. In 
our case, we do not have to worry about a secondary phonon interaction precisely because torsion does not propagate. 
On the other hand, any screening effect mediated by the graviton will require a separate, detailed assessment. 



III. COSMOLOGICAL SOLUTIONS 



In order to obtain the cosmological evolution one has to supplement the gap equation (jll.46[) and the equation 
determining the chemical potential pi. 491) with the Friedmann equation 



i^'- Y(Pgap + Pm), (III.51) 



where pm is any additional matter component. These three expressions (plus the continuity equation for pni) determine 
the evolution of the three unknown quantities A(i), /i(t), a{t). In practice, these transcendental equations are not 
analytically tractable in their full generality. However, one can solve them in some limiting regimes which can 
qualitatively capture the basic cosmological evolution. We will now focus on these regimes. In the next subsection 
we will present numerical details. 

The equation of state for the gap can be defined only implicitly via Eqs. pi.46p and (jll.50[) . We will show that there 
are different regimes where the gap scales effectively as radiation (up to a sign) or a cosmological constant, which 
confirm the expectation that the gap equation of state must be nonlinear. One can define the gap pressure by taking 
an effective Raychaudhury equation, i.e., by defining the effective barotropic index via 

1 + Wcff = , (III.52) 

6 a ma 

as a function of the scale factor. This is what dark energy experiments would measure. 



In this context we note that having a nonzero chemical potential implies that we have broken the particle-antiparticlc symmetry, which 
in turn breaks local Lorentz invariance regardless of the chosen global metric. Thus we are working on the assumption that some other 
mechanism was responsible for creating the "initial" particle-antiparticle asymmetry. 
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A. Analytical treatment 

1. Dark radiation (early times) 

To begin with, it is important to see whether the evolution of the gap and the chemical potential at early times 
can be consistent with the usual matter/radiation dominated decelerating universe. To answer this question, let us 
look into the limiting case M <^ |/^|, A and N ^ 1 (later we motivate the last condition phenomenologically) . The 
gap equation simplifies to 

A^ « 2/1^ . (III.53) 

The energy density of the gap (jll.SOP becomes 

Pgap«-32^(2Af + 3 + 21nA2) <0. (III.54) 

One can also determine the approximate behaviour of the energy density with respect to the scale factor. From 
Eq. (|II.49p . we find that /x must be negative and 

^7V+l + lnA2) = IJ , (III.55) 



or approximately (if A varies slowly) 



A~i ^ pg^p^_J_. (III.56) 

In other words, the gap energy density behaves approximately as negative radiationH This fact was already pointed 
out in 0. It is also clear that the gap energy density violates the null energy conditior0 and this property was 
exploited in @ to resolve the big bang singularity via a nonsingular bounce, as long as the equation of state of 
ordinary matter is w < 1/3. 



2. de Sitter phase (late times) 

It is relatively easy to see how a late-time de Sitter phase can emerge from the system of cosmological equations. 
If ^ A, the gap equation reduces to (/i < in order to have n> {)) 

2Tr^M^ w A^ (iV + 1 + In A^) , (III.57) 

while the total energy density is given by 



A4 

Pgap « 32^2 



(2iV + 3 + 21nA2) . (III.58) 



First, we observe that the solution to the above equation always has pgap > 0, and therefore corresponds to a de Sitter 
regime as A approaches the constant value given by Eq. (jIII.57|) . This can be seen by rewriting Eq. (IIII.58P as 



^-^"32^(,^^ + V^'- 

We are specifically interested to see whether we can explain the present dark energy driven acceleration. For this 
we require (we temporarily restore energy units) pgap ^ (meV)^ A ^ meV. When can we have such small vacuum 
expectation value for A? First of all A ~ meV corresponds to InA^ w —140. Since M"^ > 0, (|III.57p tells us 
that N > 140. Further, from (jIII.57[) it is easy to see that there are two different regimes in the parameter space 
{M,N) when we can get a small vacuum expectation value for A. If M ^ meV and 140 < N < 10^, the solution 



^ A similar contribution appears in braneworld [3C| and Hoi'ava-Lifshitz 31] cosmologies. 

® This is not very surprising, as it is well known that energy densities associated with vacuum shifts, such as Casimir energies, can indeed 
violate the energy conditions (32| . 
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corresponds to A '-^ M '-^ meV. A second possibility is to consider M <^ meV and N ~ 140. In this case we have 
A ~ e~^/^ ~ meV. In any case, the relevant range of parameters corresponds to large N and strong-coupling regime, 

> 1 , M < 1 . (III.59) 

At this point one may be concerned about the tiny value of M that is required to account for dark energy. Indeed, 
naturalness arguments would suggest M ~ 1 (Planck scale). We first point out that the M appearing in the effective 
potential (jll.45|) can be interpreted as a renormalized mass (see @ for a more detailed discussion), and therefore in 
general it can be different from the bare coupling mass in Eq. (jll.lOp . We shall not study the renormalization group 
flow of this model, and the fact that it is an effective nonrenormalizable model makes it rather difficult to interpret 
the relation between physical scales and parameters. A much deeper understanding of the quantum field theory giving 
rise to the condensate will be necessary to clarify this point, which will admittedly remain unresolved in this paper. 

We also mention that there are two other possible explanations for a small coupling M . One is tightly related with 
the fractal interpretation given below Eq. (jll.44[) . Another is to switch on a strong four-fermion interaction already 
in the Dirac action (jll.3l) . so that the gravity- induced interaction is negligible and M is actually independent from 
the Planck mass. 

Whatever the interpretation of the four-fermion interaction, and since we have not enough input to predict the 
scale at which condensation takes place, for the time being we must content ourselves to notice that a tiny value of M 
may be very compelling phenomenologically. In fact, the same gap A ~ M may also be able to account for neutrino 
oscillations which, as is well known, happen at the same mass scale as dark energy [s^ . 

The cosmological evolution of the gap energy is now clear. "Initially," A ~ ^ ^ M, and the gap energy behaves as 
negative dark radiation. Provided the very early universe is dominated by an energy density component which redshifts 
slower than radiation (weff < 1/3), such as during inflation {wcs ~ —1) or a stringy thermal phase (woff ~ 0) [34-38'], 
the negative gap energy ensures the existence of a nonsingular bounce point where the gap energy density precisely 
cancels that of ordinary matter. After the bounce, the gap energy density redshifts away faster than regular matter 
and remains subdominant as compared to ordinary matter/radiation. Thus we can have the usual decelerating phase 
of the standard cosmological model. However, once A ~ ~ M , we gradually fall into the constant gap regime 
discussed above, where A ^ A/ ^ /i. Once the matter energy density drops down to pm ^ M^, we enter the present 
dark energy dominated de Sitter phase. 



B. Numerical Explorations 



We will now verify numerically that we indeed obtain the late-time cosmology discussed above, and, in particular, 
undergo a transition from an early decelerating to a late accelerating phase. To this purpose we define 

0=-lnA2. (III.60) 

The gap equation becomes 

2t:^M^ = e'"''{N + 1 - <j>) - 2fi^{N + 2 - c/)) . (III.61) 

From Eq. (jll.49|) one obtains 



n=^ = -^e-'f'iN + l-cj,). (III.62) 



Without loss of generality we fix uq — l/(27r^). Inverting with respect to fi, 



while the gap density is given by 



Pgap = ^ (e-^ - 8m') (2A^ + 3 - 2</.) . (111.64) 



The chemical potential is negative as long a.s <f> > —1 and A^ > 0. Assuming that 4> > (the gap A decreases in time) 
ImI decreases as well, as we shall see later. 
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Plugging Eq. (jIII.63p in (jIII.61|) and solving for the scale factor, one gets 

2e3-^(iV + 2-0) 

(A^ + 1 - 0)2(iV + 27r2e'^M2) 



(III.65) 



where we chose the real positive root. By numerically inverting pil.651) we can identify two disconnected branches 
for (/)(a) as depicted in Fig. [T] 

(A) A branch confined within the interval — oo < (p < (f>A, where 0a solves the equation N +1 — (pA — 2TT'^e'^^lvP = 0; 

(B) A branch confined within the interval 0B = Af + 2<0< +oo. 

From Eq. (|II.49|) . we also observe that while in branch (A) /i < 0, in the (B) branch > 0. This suggests that while 
branch (A) is relevant for strong coupling and BCS condensation, branch (B) describes solutions in the weak-coupling 
regime where the gap is exponentially suppressed with respect to the chemical potential. 
Differentiating Eq. (jIII.65[) with respect to time, we find 



AT + 1 - - 27r2M2e<^ iV+1-0 7V + 2- 0y6 

The two branches are separated by a singularity since the Hubble rate diverges as (/> — 0a or — > (/)b- As 0a < 
+ 1 < we cannot go from one to the other. Thus, identifying the relevant branch is important for the choice 
of initial conditions, as well as the ensuing cosmology. Figure [T] shows how the branches change with respect to M 
and N: M determines 0a and the increase rate of the scale factor, while N determines only 0a and 0b. This is 
reassuring as the parameter N is arbitrary and its role is just a translation, although its actual value does determine 
the physical scale of Pgap- 

Since we are interested in the (A) branch with negative chemical potential representing the BCS condensation 
phase, we shall evolve the equations of motion from initial conditions typical of this branch. For the extra matter 
component we consider nonrelativistic dust, for the purpose of illustration. Thus we set p,„ = a~'^, M = (strong 
coupling) and = 0. These values will not correspond to the observed universe but they will capture the qualitative 
features of the dark energy solution. The acceleration of the universe is encoded in the first slow-roll parameter 

e > 1 corresponds to a decelerating universe, while e < 1 signals acceleration, if the universe expands. 

In Fig. [2] we show the evolution of the gap and its energy density in synchronous time. As — >■ 1, A — >■ e"^/2 ~ 
0.6065. The gap A decreases in time, so that the difference between the Fermi sea energy and the true vacuum of 
the theory becomes negligible. Initially, in the decelerating phase '-^ A as we expected from analytical arguments. 
However, at later times, ^ — > 0~, meaning that the density of Cooper pairs (which is positive, consistently) decreases; 
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FIG. 2. Numerical late-time solution for M — = N. From left to right and top to bottom: </>(t), gap A{t), chemical potential 
fi{t) and gap energy density pgap(t). In the last panel the matter contribution (dashed line) is shown for comparison. 



this is because the formation of pairs becomes less and less favorable. Initially pgap < and the matter contribution 
dominates. The gap density, however, increases in time and eventually (after having changed sign) dominates over 
dust. When this happens we enter the late-time asymptotic de Sitter phase. 

The evolution of the scale factor and its derivatives is plotted in Fig. [31 At early times the universe is dominated 
by the pressureless matter component pm and in fact a ^ t^/', where e ^ 3(1 + w)/2 ^ 3/2. As one can see from the 
last plot, acceleration is triggered when pgap ^ Pm but slightly before Pgap > Pm- When radiation is also added, we 
have checked that the early-time behaviour changes accordingly, e ^ 2, but the overall picture remains the same. 

The effective equation of state of the condensate, Eq. (jIII.52p . is shown in the last panel of Fig. [3l The plot starts 
from t — 0(10) because, in accordance with Fig. ^ matter dominates before that time. A characteristic prediction 
of this solution, which can be tested by future observations, is that WcS < —1 and WoS > at late times. The actual 
value of the barotropic index at the onset of the dark energy era will depend on the couplings of the model, which 
we have not tuned with respect to observations. Since we have a phantomlike regime at intermediate-to-late times, 
one should be able to distinguish our model from a pure cosmological constant or quintessence field at high enough 
redshift. The fit with the supernova data goes beyond the goals of the present investigation. 



IV. DISCUSSION 

In this paper we have proposed a mechanism realizing late-time cosmic acceleration due to the possibility that 
fermions can condense in the very early universe. We have shown that all of the necessary conditions for condensation 
exist, including a covariant attractive four-fermion interaction. By analyzing the transcendental equation relating 
the gap and the scale factor we discovered that at late times the contribution of the gap energy drives a phase of 
acceleration. The details of the model depend on the regularization scheme and the value of the couplings, but the 
emerging qualitative picture is robust. 

In this context, we note that one may be worried that, given the low value of M required for phenomenology, 
M meV, our procedure for integrating out gravitational degrees of freedom is not justified, and whether we might 
be missing corrections of order A/M, which may be important as A/M is not small. Let us reemphasize why the 
method we adopted to obtain the effective potential for A is robust. In the first step, we integrated out the torsion 
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FIG. 3. Numerical late-time solution for M = — N. From left to right and top to bottom: a{t), Hubble parameter H(t), first 
slow-roll parameter e{t) and effective barotropic index Wcsit) of the condensate equation of state. 



field to obtain a four-fcrmion interaction. This was done by imposing the classical equations of motion for the 
torsion field. However, we note crucially that the torsion is a nondynamical field: its field equation is algebraic and 
there are no derivatives acting on the torsion. In other words, the torsion field does not propagate. If one were to 
draw Feynman diagrams involving torsion propagators, these would just collapse into points leading to four-fermion 
interaction vertices. Unlike gauge theories where the four-fermion interaction is only an "effective" interaction valid 
below the mass of the gauge bosons, the torsion-mediated four-fermion interaction is an exact quantum mechanical 
description, because the torsion field has no kinetic term. Therefore the actions given by (jll.ip and (jll.3|) are exactly 
equivalent to (jlLlOp (plus the usual Einstein-Hilbert and Dirac actions) both classically and quantum mechanically. 
This equivalence does not depend on how small or large the value of M is. 

In the next step, we integrated out the fundamental fermionic degrees of freedom to obtain an effective theory of 
their bound state, the gap A (Cooper pairs). The procedure we followed is precisely the same that is adopted in 
condensed matter literature. In particular, this is known to capture not only the weak-coupling BCS limit {M — cx)), 
but also the strong-coupling BEC limit (Bose-Einstein Condensate, M — )■ 0), which is what we are interested in 
(for general reviews on BEC and the BEC-BCS crossover, see |39l - t43| ). We note that the effective potential for A 
features a ln(A/Af ) term, which contains all the powers of (A/M). In this sense, the effective action Scs in Eq. (jll.33p 
describes nonperturbative phenomena by summing all such terms. 

Having clarified the nonperturbative character of the condensation, it remains to explain the small value of M, 
which is a partially unrelated question. We have not answered it here but we have commented on several possibilities, 
including a renormalization mechanism and the presence ah initio of a strong fermionic interaction. To summarize, one 
is indeed assuming that a bare C(l) cosmological constant "miraculously cancels" by virtue of some mechanism. So, 
we do not solve the "old cosmological constant problem" and in this sense fair no better than the usual quintessence 
models, where one assumes that all contributions to the effective cosmological constant conspire to reduce it to 
a dynamical field with a "small" potential. However, our model has some distinct advantages. First, because the 
suppression of the cosmological constant is exponential, the level of fine-tuning is reduced to just one part over 100, via 
the choice of N . Therefore, if one could give a physical interpretation to N , then the "smallness problem" associated 
with the observed value of dark energy would be relaxed. This was attempted in Sec. IIIDI Second, our situation is 
consistent because we have calculated the nonperturbative potential, whereas in usual quintessence models it is not 
clear why higher-order corrections are suppressed. 
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There are several issues we have not considered. Since we have worked on the assumption that, hke in Minkowski, 
the fermion interaction is nonrenormahzable, the regularization parameter N of the effective theory has been assumed 
to be physical and therefore should ultimately be motivated by the fundamental microscopic theory. Also, we have 
assumed to be a constant. A possibility we leave for future study is to allow for a time-varying N{t) or coupling 
M{t). For instance, previous literature [HI, [l^j El] have considered different choices for a physical cutoff, X — fi or 
X — H , which would correspond to having a time-varying aQ This may lead to interesting situations, including the 
possibility of not having to tune M to the very tiny meV scale. 

In order to verify the robustness and observational validity of the rich cosmological picture we presented here, 
future studies will have to go into greater detail in the analysis of the parameter space and initial conditions. It is 
promising that the same condensation mechanism can solve the big bang singularity and the dark energy problem. 
We end by pointing out that the dark energy scenario presented above has some distinctive observational features. 
For instance, depending upon the detailed history of the early universe the negative dark radiationlikc gap energy 
may be detectable in BBN and CMB observations (sil. Is^j. Also, there is a most encouraging possibility of linking 
the scenario with neutrino physics d, as mentioned^a^iTer. 
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Appendix A: Derivation of the four-fermion interaction 



The Hoist action can be written as (G = 1) 



where 



= ^ / d^xee^e'^jP'^LRj^ , (A.l) 



P"kl = S{i,Si^ - ^^"kl, (A.2) 



-1 KL _ 7 xK xL , ^ ^ KL 



whose inverse is 

P-\j^^^-^^[8f,8':, + -e,,-^y (A.3) 
Variation of the Hoist action with respect to the connection yields 

Likewise, variation of the Dirac action gives 

5Sd 



l-D,(ee\^e''])p"KL. (A.4) 



5AKL 8 



eV'{7[-R'7L],7^}eyV' 



^e'A'LM(^757*V)e^ (A.5) 



In the second line we have used the identity {'^[kIl] ; 7/} — '^^^iKLMlbl^'^ ■ The total variation of the action 5h + 
with respect to the connection is 

(e el^e;;!) P'' kl = 2tt et' klmJ^' e'i , (A.6) 



The ambig uity in the choice of the cutoff is not dissimilar to the one entailed in modern and inequivalent formulations of loop quantum 
cosmology |45|, l46l . There, one can choose the kinematical area of the elementary holonomy to be fixed in time and equal to the Planck 
area (improved quantization scheme 53)' °^ ^'^^ make it dynamical as in lattice refinement models [4St45Clll , which corresponds to 
probing geometry with a time-dependent microscope. 
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where J^^ = 'tp'ys'y^ip is the axial current. 

Writing the connection as Ajf = ujj/ + Cj/ , where oj is the connection compatible with the tetrad, and using 
Eq. (|A.6p . one gets 



^^MP'^Q] + C^pq: = 2vr-^er J5M U" pq + ijf/^^, j . (A.7) 



7^ + i V 7 

Contracting with we obtain 



From the above two equations we obtain 

v,2 



72 + 1 



(^MK"ef-^^i.^,ej). (A.9) 



Inserting the above expression into the first-order gravity + matter action yields the four-fermion interaction. For 
the calculation in the Nieh-Yan case we refer the reader to [13, ■ 
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